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Passive Suppression of Pogo on the Space Shuttle 

by 


M. H. Lock and S. Rubin 
The Aerospace Corporation 


ABSTRACT 


This study provides a qualitative assessment of (1) the tendency for 
pogo instability of the Shuttle vehicle in the absence of suppression devices 
and (2) the effectiveness of two passive suppressors (the compliant accumu- 
lator and the resistive accumulator) in counteracting any tendency toward 
Instability. In addition, the relative effectiveness of three suppressor loca- 
tions [the low pressure pump (LWOP) inlet and discharge and the high pres- 
sure pump (HPOP) inlet] is also evaluated. 

The primary conclusion of the study is that effective pogo suppression, 
with passive devices, can be accomplished at the HPOP inlet location but 
not at the LPOP locations. Other conclusions are drawn regarding the 
relative effectiveness of the two accumulator types and with respect to tenta- 
tive accumulator design requirements. Finally, a number of recommenda- 
tions are made regarding future studies. 
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PASSIVE SUPPRESSION OF POGO ON THE 
SPACE SHUTTLE 

by M. H. Lock and S. Rabin 
The Aerospace Corporation 


SUMMARY 


This study provides a qualitative assessment of (1) the tendency for 
pogo instability of the Shuttle vehicle in the absence of suppression devices 
and (2) the effectiveness of two passive suppressors (the compliant accumu- 
lator and the resistive accumulator) in counteracting any tendency toward 
instability* In addition, the relative effectiveness of three suppressor loca- 
tions [the low pressure pump (LPOP) inlet and discharge and the high pres- 
sure pump (HPOP) inlet] is also evaluated. 

The assessments are made upon the basis of stability analyses under- 
taken with an idealized model of the Space Shuttle. Although an idealization, 
this model contained description of the higher organ-pipe modes of the feed- 
line, an inter pump mode of the oxidizer system, and the coupled longitudinal 
and lateral motions of the vehicle structure. The model is especially rigor** 
ous with regard to the feedline modes in that a continuous representation is 
employed. In addition, the generalized force contributions due to fluid 
resistance, fluid -momentum effects , turns in the fluid flow, and tank out- 
flow effects are considered as well as the engine thrust contributions. 

The results of the stability analyses revealed instabilities of the basic 
system (i.e., no suppressors) that involved the inter pump mode as well as 
the feedline-type modes of the propulsion system. Both instabilitie s were 
found to be eliminated by the use of either compliant or resistive accumula- 
tors at the HPOP inlet. The use of these accumulators at either the LPOP 
inlet or discharge was found to aggravate the instability that involved the 
Lnterpump mode. Accumulators at the LPOP inlet were found to eliminate 
the instabilities that involved the feedline-type modes of the propulsion 
system. 

The primary conclusion of the study is that effective pogo suppression, 
with passive devices, can be accomplished at the HPOP inlet location but not 
at the LPOP locations. Other conclusions are drawn regarding the relative 
effectiveness of the two accumulator types and with respect to tentative 
accumulator design requirements. Finally, a number of recommendations 
are made regarding future studies. 
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INTRODUCTION 


The suppression of pogo instability on the Space Shuttle is complicated 
by the nature of the vehicle structure and liquid propulsion system. In the 
case of the structure, the multibody configuration introduces significant 
coupling between the lateral and longitudinal motions of the vehicle. Since 
all such motions can couple with the propellant flows, the analysis of the 
pogo phenomenon on the Space Shuttle will require an accurate description of 
both lateral and longitudinal motions. In the case of the propulsion system, 
the length of the oxidizer feedline results in a number of higher organ-pipe 
hydraulic modes in the frequency range where strong structural/propulsion 
system interaction could occur, thereby raising the possibility of pogo insta- 
bility over a broad range of modes and frequencies. In addition, the engine 
itself can introduce a propulsion-system mode that is primarily internal to 
the engine and which results from the presence of two pumps with a signifi- 
cant length of intermediate ducting. An unpublished study by S. Rubin has 
suggested that such an "Lnterpump" mode of the oxidizer system can lead to 
instability. This latter circumstance raises a question about the location of 
a pogo suppression device since it is possible that a suppressor located at 
the inlet of the engine would not be effective in suppressing an instability that 
involves an Lnterpump mode. Conversely, a suppressor located internal to 
the engine may not be effective in suppressing an instability that involves a 
feedline mode. 

The goal of the present study is to provide judgments on the issues 
raised by the complexities of the Space Shuttle system. These judgments 
are developed from an analysis of the system which, although simplified, 
does contain an account of the lateral and longitudinal motions of the struc- 
ture, the higher oxidizer feedline organ -pipe modes, and an Lnterpump mode. 
Specifically, the study is directed at (1) an examination of the tendency for 
pogo instability of the Shuttle vehicle in the absence of suppression devices, 
(2) an evaluation of two passive accumulator types, compliant and resistive, 
and (3) an evaluation of three accumulator locations. 

The compliant accumulator (no inertance or resistance) is a simplifi- 
cation of the type of device commonly used on past vehicles. These devices 
were, in general, intended to be purely compliant but, in fact, possessed 
some small inherent inertance and resistance. The resistive accumulator 
(substantial resistance, along with inertance and compliance) was introduced 
quite recently on the Delta Stage I vehicle to help correct a problem that had 
arisen due to significant spatial separation between the accumulator and the 
engine inlet (ref. 1). This separation introduced a second hydraulic mode of 
the feed-system that led to instability. 

The consequences of the Delta Stage I instability were ameliorated with 
the introduction of significant resistance by means of the accumulator. 


A companion study (NAS3- 17758) is being undertaken on the use of active 
suppressors. 
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Although the device had been introduced for a problem that resulted from 
accumulator /pump separation, the character of the device indicated that it 
could provide system damping over a wide frequency range irrespective of 
accumulator /pump separation. This feature of the resistive accumulator 
made it appear an attractive candidate for use on the Space Shuttle and led to 
its selection for examination in the present study. For simplicity, the 
present study only treated the case of zero accumulator /pump separation. 
The examination of the ' conse quence s of finite separation (applicable for both 
resistive and compliant devices) is considered more appropriate for later 
more refined studies. 

With respect to the evaluation of accumulator location, it was decided 
to perform stability analyses with both accumulator types at the following 
three locations: 

1. Low-pressure pump inlet 

2. Low-pressure pump discharge 

3. High-pressure pump inlet 

The results obtained for these cases were also compared against the cor- 
responding results for the basic system (i.e., the vehicle without suppres- 
sors). The analyses were undertaken for the oxidizer system alone since 
unpublished preliminary studies by the Rockwell International Space Division 
have indicated that this portion of the propulsion system is most likely to be 
critical from the standpoint of pogo instability. 

The subsequent sections of the report describe the development of 
design criteria for the resistive accumulator, the formulation of the ana- 
lytical model for the stability analyses, and the stability analyses of the 
basic system and the system with incorporation of the compliant and resis- 
tive devices. The results are summarized, conclusions are drawn, and 
recommendations made for future studies. 

The authors would like to acknowledge the work of Raymond E, Orth 
and Heather Bagwell in the programming of the system stability equations 
and the propulsion system equations, respectively. 
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NOMENCLATURE 


a 

A 

A, 

B 

C 

C f 

D 

T 

g 

G 


L 

m+1 

M 


M 


n 




n 


n 


Q 


Units; M (mass), F (force), L (length), T (time) 

acoustic velocity [ LT 
area [L^] 

effective area of thrust chamber [L^] 

see eq. (1. 3) 

compliance, [F *L^] 

thrust coefficient, [— ] 

see eq. (1. 3) 

mean mass flow of propellant, [MT' 1 ] 

_2 

standard acceleration of gravity, [LT" ] 
structural gain for engine motion, 0 (e)/M, [M ] 
height of propellant in tank, [L] 
height of tank -to -engine feedline, [L] 
imaginary unit, n/^T 

line length, [L] 

r -5 2 -4 

mertance, [FL T or ML ] 

pump dynamic gain, [— ] 

mass, [M] 

generalized mass of nth structural mode, [M] 

modal tank -bottom pressure per unit acceleration of 
generalized coordinate, [FL“ 3 T 2 ] 

-2 

oscillatory pressure, [FL-‘J 

generalized displacement of nth structural mode 

3 1 

volumetric flow, [L T ] 
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t 

T 

x 

25 

V 

y 

z 


a.. 


L 

% 

e 

p 

r 

0 


’N 


n 




w 


resistance, [FL’^T] 

Laplace variable used to denote the complex frequency 
cr + iw, [T - 1] 

time, [T] 

thrust, [F] 

structural displacement along longitudinal axis, [L] 
structural displacement along lateral axis, [L] 
steady flow velocity of propellant [LT"\l 
flow admittance, [F 
flow impedance, [FL ^T] 

coefficients in feedline transmission function; see 
Appendix A [— ) 

ratio of critical damping for structural mode, [—] 

ratio of critical damping for accumulator [— ] 

structural damping required for neutral system stability, [— ] 

complex propagation angle 

propellant mass density, [ML ] 

travel time in a hydraulic line, /LC or //a, [T] 

modal displacement, [— ] 

angular frequency, [T *] 

natural frequency of the accumulator [T~^] 

-1 


u> n natural frequency of structural mode, [T~ ] 

ft frequency ratio; see eq. ( 1 . 5 ) [— ] 

Subscr ipts 

a accumulator 

b bubble 

c chambe r 
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d 

e 

ip 

l 

n 

P 

R 

t 

u 


discharge 
engine 
inter pump 
lower limit 

nth structural -system mode 
pump 

real part; relative 
tank 

upper limit 
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1. DESIGN CRITERIA FOR RESISTIVE ACCUMULATOR 


The objective of the use of the resistive accumulator is to produce a 
significant level of hydraulic damping over a specified frequency range. 

This section of the report provides a measure of the required level of damp- 
ing and develops accumulator design criteria for achieving the desired 
characteristics. 


1. 1 Accumulator Model 


The accumulator is assumed to be a linear device characterized by an 
inertance, L a , a compliance, C a , and a resistance, R a . Assuming time- 
dependent perturbations of the form e st , where s is the Laplace variable, 
the equations governing the accumulator can be written 

sLQ +RQ = P -P ' (la) 

a a a a a 

Q = - s C P o (lb) 

a a a 

where Q is the volume flow out of the accumulator, P is the accumulator 
internal pressure, and P is the pressure at the entrance to the accumulator 
(see sketch). 



The flow admittance, 



of the accumulator can be written 


Y 

a 




2 

+ 2 £, u s + w 
a a a 


-1 


(2) 


where the accumulator natural frequency, to , and damping ratio, 
defined by a 


to 


— i • r 

L_ C„ ' a 


R C to 
a a a 


V are 
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1,2 Approximate Stability Analysis 


The basis for the design criteria for the resistive accumulator was 
developed from an approximate stability analysis of the single -pump system 
shown in figure 1. The analysis was undertaken under the assumption that 
the structural system responds in a single mode. Employing the approximate 
method of reference 2, the structural damping required for neutral stability, 
is found to be given by 




(m-J-T )BG 


R Real Y +(m+l) 
e a 


(i-d; 


where 


i \ 

A 1 A T R c l 1 + h 0 

V S 


A l R e *e 

(m+l)A T R c 


and where 

A^ = cross -sectional area of feedline 

h = length of feedline 

s 

hj. = height of fluid in tank 

R = engine resistance 

c 

R = thrust -chamber resistance 

c 

Aj = effective area of thrust chamber (i.e., area x thrust 

coefficient) 

(m + 1) = gain of pump 

= modal displacement of engine 

#£ = modal displacement of center of gravity of propellant in tank 

= frequency of the structural mode involved 
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PUMP (gainm+1) 


Figure 1. Model for accumulator design criteria analy 




The structural gain, G , for longitudinal motion of the engine is defined by 



where M is the generalized mass of the structural mode. 


1.3 Design Criteria 


The form of equation (3) indicates that the structural damping required 
for neutral stability will decrease with increase of the real part of the 
accumulator admittance. The equation also indicates that the magnitude of 
this term should be such that 

Real Y >>(m+l)/R (4) 

a e 


for the resistive accumulator to be effective. Therefore, the design goal for 
the accumulator is to maintain the inequality [eq. (4)] over the frequency 
range of interest. 


For the case of harmonic oscillations at the frequency, co, the real part 

of the admittance, Y , can be written 

a 


Real Y 

a 


w L 


a a 



+ ( 2; a ) 2 


-1 


(5) 


where the frequency ratio £2 = w/u . The derivative of Real Y a with respect 
to £2 is a 


Y a = 


4C d-n 4 ) 


CO 


L £2 3 (l/£2-£2) 2 + (2£ 2 ) 
a a 


(6) 


The results [eqs. (5) and (6)] indicate that Real Y a increases monotonically 
to a maximum at £2 equal to unity and then decreases monotonically as £2 
increases above unity. The maintainance of Real Y a above some selected 
value, say o', over the frequency range, < co u> a , is assured if 

Real Y (w- ) > <* (7a) 

Sl 

Real Y (u> ) > a (7b) 

a' u 
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TVio rr\r\^ifirmc fan f 7 \] ^hsn Is to the follow in g i* c ui r c one n t c 
accumulator parameters 


u> . C yfl a ( 1 - R a) 2 
l a a a 


ui C J R a (1 - R.u ) 2 
u a a a 


oil - u > 2 L C ( 8 a) 

‘ i a a 1 

u |l - w 2 L C ( 8 b) 

i u a a I 


The requirements [eq. ( 8 )] constitute the design criteria for the re- 
sistive accumulator. For a specified design goal these inequalities define a 
region of the accumulator parameters within which the design goal is satis- 
fied. The application of the criteria is illustrated below for the design goal 
that was used subsequently in the study. This goal was 

Real Y > 0. 06 m 5 /MN s (1 in. 2 /se c ) for 2 < £ < 30 Hz * (9) 

a 


The allowable ranges of the accumulator parameters are most readily 
obtained by specifying the accumulator compliance, C , and then using 
equations ( 8 a) and ( 8 b) to define an allowable region o^inertance , L a , and 
resistance, R a . A representative range of the accumulator compliance was 
taken to be from 0. 006 m^/MN (0. 1 in. 2) to 0. 12 m^/MN (2. 0 in.' 2 ). This 
range was determined by enveloping compliance values that were estimated 
for a representative 0. 057 m3 (2 ft^) volume accumulator using propulsion 
system pressure data supplied by Rockwell. The resulting regions of allow- 
able values of inertance and resistance are shown in figure 2 together with a 
smaller region that was calculated for the condition Real Y 2 0. 3 m^/MN s 
(5 in. 2 /sec), also for 2 if i 30 Hz. From the figure, it is a seen that the zones 
are relatively insensitive to the value of the accumulator compliance (at 
least for the range that was considered) but are quite sensitive to the selected 
level of the admittance. Insofar as an accumulator design is concerned, the 
selection of any combination of inertance and resistance values that are 
within the zone defined by curve A (fig. 2) will lead to satisfaction of the 
design goal [eq. (9)]. The particular combination that was selected for use 
in the stability analyses was 

R = 1. 7 MN s/m 5 (0. 1 sec/in. 2 ) 

ct 


L = 0. 023 MN s 2 /m 5 (0. 00133 sec 2 /in. 2 

3 . 


These values were judged to be practically achievable on the basis of Delta 
Stage I vehicle experience (ref. 1). In addition to this particular design, a 
number of other resistive accumulator designs were briefly examined to 
determine the sensitivity of the stability results to the accumulator para- 
meters . 

The admittance in SI units is based upon volume flow divided by pressure; 
in English units, weight flow Is employed. This flow difference also applies 
to the resistance and inertance units. 
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2. ANALYTICAL MODEL FOR STABILITY ANALYSIS 


2. 1 Shuttle Propulsion System /Structural Model 


For the purpose of the stability analysis, the most important features 
of the vehicle propulsion system were represented using a single equivalent 
engine. The elements employed are shown in figure 3 and consist of a lox 
tank, two feedline segments (one longitudinal and one lateral), a low-pressure 
pump (LPOP), an interpump line, a high-pressure pump (HPOP), a discharge 
line, an injector, and a combustion chamber. For the purposes of the 
present study, the motion of the Shuttle vehicle was represented by a single 
structural mode. This representation was considered adequate except for 
the exceptional case of proximate structural modes (i.e., modes with very 
close frequencies) with comparable gains. As a consequence of this repre- 
sentation, the structural modes selected for use in the stability analysis 
were examined one at a time. 

The detailed equations that govern the motion of the system are pro- 
vided in Appendix A of this report. In order that these equations provide a 
representative description of the higher organ-pipe modes of the feedline, 
the exact solutions for a continuous hydraulic representation (including re- 
sistance) were employed to develop feedline transmission functions. This 
results in the appearance of hyperbolic functions of complex argument in 
the system equations. Since the wave transit times associated with the inter - 
pump and discharge lines were relatively short compared to the structural 
response times of interest, the fluids in these latter lines were treated as 
incompressible. Generalized force contributions that result from fluid re- 
sistance as well as from the convective derivative (v3v/9x, where v is the 
velocity of the flow) that appears in the fluid momentum equation were also 
included in the system equations . 

On the assumption that the structural modes are developed with closed- 
bottom tanks , tank outflow effects are included (ref. 3). All corner forces 
at turns in the flow and drag forces across resistive elements are also in- 
cluded. 


2, 2 Computational Procedure 


For the purpose of the stability calculations the equations of motion 
were reduced to a 14 th order system with the form 


[V(s)] + [E][F(s}]j H = 0 


( 10 ) 


where the matrix [V(s)] describes the basic coupled structure /propuls Lon 
system (i.e., the system in the absence of accumulators); the matrix [E] 


13 




Figure 3. Lox system elements of system stability model 
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f ir o vide s the specification of the accumulator location while the matrix 
F(s)] contains the description of the accumulator characteristics. The 
solution vector H comprises 


A 

H 


P 2' P 4’ P 5* P 7' 


P 8’ P < 


Q t ,Q 2 ,Q 3 ,Q 4 ,Q 5 ,Q 7 ,Q 8 ,q 


m 


(ID 


where P- and Q; denote the pressures and flows at various points within the 
oxidizer system (fig. 3). The q n term is the generalized coordinate asso- 
ciated with the nth structural mode; the motion 1c of the vehicle at some point 
r is related to q n by the following: 

x = q e St 0 (r) 
n N 


where $ n fr) is the structural mode shape and s is the Laplace variable. 

The matrix equations are completely defined when the propulsion- 
system parameters and structural -mode data are specified. A When this is 
done, the eigenvalues, s^, and corresponding eigenvectors, H^, of the sys- 
tem can be obtained by satisfaction of the conditions 

det |[V(s.)] + [E][F(s.)]|= 0 (12a) 

j[V(s.)] + [E][F(s.)]| H. = 0 (12b) 


Because of the use of the exact feedline solutions the determinantal 
equation is a transcendental equation in the Laplace variable. The eigen- 
values of the equation were obtained with the use of an iterative root-finding 
subroutine (ref. 4) that used the input structural frequency, w n , and pre- 
viously calculated propulsion-system eigenvalues as initial guesses. 


2. 3 Numerical Values of System Parameters 


2. 3, 1 Propulsion-System Parameters 

The cross-sectional areas and lengths of the various lines were 
based upon available Space Shuttle design data. The resistance and inertance 
of the lines, pumps, engine, and thrust chamber were developed from the 
basic data given in the SSME Engine Dynamic Model (ref. 5) and from infor- 
mation in a previous pogo study (ref. 3). The estimated values of these 
parameters are provided in Appendix B. The variation, with time of flight, 
of the cavitation bubble compliance at the pump inlets and of the pump gains 
were estimated from available operating data; the compliance estimates were 
made using the results presented In reference 3 (obtained from the ’’stay- 
time" method of ref. 6) and the pump gain from unpublished Titan and Delta 
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vehicle studies. The resulting time variations are shown In figures 4 and 5. 
The compliance and gain values employed in the stability analyses at speci- 
fic flight times were taken from these figures. 

2.3.2 Structural-Mode Data 

The structural-mode data employed in the stability analys es‘ were 
taken from a set of modal calculations undertaken by Rockwell. These cal- 
culations had been performed for the case of symmetric vibrations of the 
vehicle (i.e., no yaw motion) and were based upon a vehicle configuration 
(designated M89B) that was current in early 1973. The data were considered 
to be the best that were available for the present study. It may be noted in 
passing, that the Space Shuttle vehicle has subsequently been reduced in size 
and weight; however, modal data for the new configuration were not available 
for this study. The Rockwell data comprised the frequencies and mode 
shapes of the first hundred symmetric vibration modes at five specific flight 
conditions; in addition to the structural data, the lox tank-bottom displace- 
ments and pressures were also provided. The set of conditions treated by 
Rockwell is given in the table below together with the associated times of 
flight and the corresponding maximum and minimum values of the calculated 
frequencies; the abbreviation SRB that appears in the table denotes the solid - 
rocket booster. 


Condition 

Flight Time 
(sec) 

Frequency (Hz) 

h 

100 

Liftoff 

0 

2. 18 

44. 5 

Max. Dynamic 
Pressure 

54 

2.24 

48.9 

Before SRB 
Separation 

116" 

2.28 

53.3 

After SRB 
Separation 

116 + 

2. 29 

62. 4 

Orbiter 

End-Burn 

480 

2. 81 

103. 8 


Examination of the modal data from the standpoint of structural 
gain indicated that the higher gain levels were associated with modal fre - 
quencies above 20 Hz. Based upon this feature of the calculated results 
and the increasing uncertainty of the data for the higher modes, a frequency 
of 3 0 Hz was selected as a reasonable upper limit to the frequency range to 
be considered in the subsequent stability analyses. The data also indicated 
that the higher levels of the structural gain, in the frequency fange of inter- 
est, tended to be relatively insensitive to time of flight. This feature is 
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CAVITATION COMPLIANCE, 10" J rrr/MN 
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0. 1r* 


1 


-L 
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-U 
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0. 001 


Figure 4. Variation of pump inlet cavitation compliance 






Figure 5 



ariation of pump gain 




illustrated in figure 6 where the variation with flight time of the two 
largest structural gains is presented; for comparison of relative magni- 
tudes the gain in the fundamental mode is also shown. In the figure the 
higher gains have been identified by the approximate frequency of the struc- 
tural mode. With the exception of one case, the gains are associated with 
a single mode. The exceptional case occurred at the orbiter end-burn 
condition where two proximate modes, one at 26. 8 Hz and the other at 
2 7.2 Hz both exhibited large structural gains of comparable magnitude. 


For application in the stability analysis, the calculated structural- 
mode frequencies were allowed to vary through some ±15 percent. This 
variation was introduced to provide at least some account for changes in 
vehicle configuration and also uncertainties in the structure /propuls ion- 
system modeling. The variation was expected to cover the worst case 
conditions in terms of the proximity of structural and propulsion resonances. 
To account for damping in the vehicle, a critical damping ratio of 0. 01 
was assigned to each structural mode. 


Since the stability analysis is based upon a single equivalent system 
with three identical engines, the engine modal amplitude, 9 S e> used in the 
stability calculations is related to the modal data by 

3 



i= 1 


where the are the modal amplitudes of the individual engines. In a simi- 
lar manner the tank-bottom pressure excitation.^ used in the analysis is 
related to the corresponding modal data by 


& 

n 




3. SYSTEM STABILITY ANALYSIS 


To initiate the study, the following three suppressor configurations 
were analyzed: 

1. Basic system (i.e., no pogo suppressor) 

3 3 

2. System with 0. 057 m {2 ft ) volume compliant accumulator 

3 3 

3. System with 0. 057 m {2 ft ) volume resistive accumulator 
designed to meet the requirements: Real Y a 5 0, 06 m^/MN s 
(1 in. ^ /sec) for 2 < f < 30 Hz. 
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The 0. 057 m 3 (2 ft 3 ) volume used for both the resistive and com- 
pliant accumulators was chosen as a representative practical value. The 
condition on Real Y a for the resistive accumulator was based upon an esti- 
mated magnitude of 0. 008 m^/MN s {0. 14 in. 2 /sec) for the comparable pro- 
pulsion-system term (m + l)/R e . Thus, the design ensures that Real Y a >> 
(m+l)/R_ over the frequency range of interest. 


3. 1 Basic System 


3. 1. 1 Propuls ion -System Modes 

The eigenvalues and eigenvectors of the isolated propulsion sys- 
tem (Le., stationary structure) were calculated both to provide initial 
guesses for the root-finding program used in the subsequent stability analysis 
and to identify the interpump mode of the propulsion system. The calcula- 
tions were made using the numerical data given in Appendix B together with 
cavitation compliance and pump gain values taken from figures 4 and 5; the 
values of these latter items for the orbiter end -burn, liftoff and after SRB 
separation conditions are given in the following table. 


Case 

Item 

Liftoff 

After SRB 
Separation 

End-Burn 

LPOP Gain 
(nrij+l ) 

1. 625 

2. 2 

1. 306 

HPOP Gain 

( m 2 +l ) 

1.48 

1. 54 

1.42 

LPOP Inlet, C L1 
b 1 

(10' 3 m 5 /MN) 

1. 2 

2. 3 

0. 59 

(in. 2 ) 

(0. 0204) 

(0.039) 

(0. 01) 

HPOP Inlet, C b2 
(10- 3 m 5 /MN) 

0. 27 

0. 3 

0. 24 

(in. 2 ) 

(0. 0045) 

(0. 0051) 

(0. 004) 


The modal frequencies and the associated critical damping ratios calculated 
for the orbiter end-burn condition are presented graphically in figure 7. 

The results for the other flight conditions are similar. The interpump mode 
is noted in this figure. It is seen to be relatively highly damped with a 
frequency of 24.4 Hz (the frequency of this mode at the liftoff and after SRB 
separation conditions was 23. 1 and 21.6 Hz, respectively). This resonance 
can be viewed simply as involving the interpump and LPOP fluids moving on 
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the LPOP inlet cavitation bubble compliance at the upstream end and the 
HPOP inlet cavitation bubble compliance at the downstream end (fig. 3). 
Consequently, the resonant frequency can be estimated with the formula 


£. (Hz) Rs 

Lp 


2tt 



bl 


+ 1/C 


L, 4- L . 
L P l 


b2 


1/2 


where L^, Lp^ denote the inertance of the interpump line and the LPOP, 
respectively; and where C b 1 ' C b2 denote the cavitation bubble compliance 
values at the inlet to the LPOP and HPOP, respectively. The above expres- 
sion provides a means of gaging the sensitivity of this modal frequency to 
changes in the involved parameters. The remaining oxidizer system modes 
can be considered as “feedline -type ” modes since their frequencies lie within 
the frequency bands defined by the open-open and open-closed modes of the 
feedline. The different character of the modes is illustrated in figure 8 
where representative fluid pressure amplitude distributions are shown for 
the interpump mode and the first two feedline -type modes; the in or out 
phasing, shown crudely, approximates the actual phasing. The appearance 
of significant amplitudes only within the engine (i.e., beginning with the 
LPOP) for the interpump mode is clearly seen from the figure. 

3. 1.2 Selection of Cases 


A particular stability case for the basic system comprised a speci- 
fied flight condition and a specified structural mode. For a given case, the 
stability calculations were performed over the assumed ±15 percent varia- 
tion in the structural mode frequency. Selection of the stability cases, 
particularly the selected structural modes, was largely governed by the 
character of the modal data; the modes with the higher structural gains 
being the candidates for the analysis. The end-burn condition was selected 
for initial analysis. The set of modes employed for analysis of end-burn 
stability are given in the table below, with the associated modal frequencies 
and the structural gains. In the table the nomenclature E denotes end-burn 
while the numbers give the order of the mode. 


Mode 

Item 

El 

E2 

E7 

E30 

E34 

E35 

f n (Hz) 

2. 8 

4. 7 

8. 5 

22. 5 

26. 8 

27. 2 

g (itr 6 /k g ) 

tJ 

1. 25 

1. 14 

6. 3 

57 

177 

274 

(10’ 4 in. /lb-sec^) 

(2. 2) 

(2. 0) 

(11) 

■■ —J 

(100) 

(310) 

(480) 
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Figure 8. Lox system modes: orb iter end-b 




It should be noted that the structural gains assigned to the thirty -fourth and 
thirty-fifth modes are twice the values given by the modal data. This was 
done to account for the possibility of destabilizing coupling between the modes. 
Such coupling could not be described by the present analysis since it involved 
only a single structural mode. The use of a factor of two did not seem un- 
reasonable since the two structural gains are of a comparable order of mag- 
nitude . 


The tank-pressure excitation at the orbiter end-burn condition 
was assumed to be negligible. Since the absence of this excitation at end-burn 
could tend to de -emphasize the role of the feedline -type modes, it was clear- 
ly necessary to treat at least one earlier flight time where this excitation was 
present. The condition initially selected was liftoff where six structural 
modes were considered in the stability analyses. These modes with their 
associated frequencies and structural gains are given in the table below 
together with the corresponding values of the modal tank-bottom pressure, 
per unit engine acceleration, /$ e * 


Mode 

Item 

LI 

LI 8 

L26 

L49 

L51 

L62 


2.2 

8.9 

12. 1 

21.9 

22.5 

27. 1 

G (10 _6 /kg) 

9. 1 

4.6 

11.4 

74 

8.6 

206 

(10 4 in. /lb-sec 2 ) 

(16) 

(8) 

(20) 

(130) 

(15) 

(360) 

(Ns 2 /m 3 ) 

115 

5.7 

3. 0 

820 

7640 

-1.0 

(10 ^ lb -sec^ /in. ^ ) 

(42) 

(2. 1) 

(1. 1) 

(300) 

(2800) 

(-0.4) 


Subsequent to the performance of the analyses at liftoff, two structural modes 
for the after SRB separation event were also selected for use in stability 
analyses at this latter flight time. These modes and their characteristics 
are given below. 
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Mode 

Item 

A 1 

A46 


2. 3 

27. 0 

G e (10" 6 /kg) 

7. 4 

194 

(10" 4 in. /lb-sec 2 ) 

(13) 

(340) 

/<f> (Ns^/m 3 ) 

n e 

464 

-0. 014 

(10~^ lb-sec^ /in." 3 ) 

(-170) 

(-0. 005) 


3. 1. 3 Stability Results 

3. 1.3.1 Orbiter End-Surn 

The stability results obtained for the basic system at the orbiter 
end-burn conditions are presented in figure 9. The calculated system criti- 
cal damping ratios are shown as functions of the assumed structural natural 
frequency when toleranced ±15 percent. Recalling that the damping of the 
structural mode is input at 1 percent of critical, it is seen that the damping 
in the first and second structural modes (f ~ 3 Hz and 5 Hz, respectively) is 
relatively unaffected by the action of the propulsion system. The seventh 
and thirtieth modes (f ~ 8 Hz and 22 Hz, respectively) show a somewhat 
greater effect; however, the minimum value of the calculated damping ratios 
remains above 0. 007. The thirty-fourth mode (f ~ 27 Hz) is strongly affected 
with the minimum damping ratio being less than 0. 002 while the thirty -fifth 
mode (f ~27 Hz) is the worst case with instability being predicted if the 
structural mode frequency lies within the range 23 to 25 Hz. The damping 
ratio in this range reached -0. 001. In this last case, the propulsion- 
system variables in the unstable mode (fig, 10) contain contributions from 
the inter pump mode (real part of the mode shape) and from the sixth feedline- 
type mode (imaginary part of the mode shape); however, examination of the 
energy transfer from the propulsion system to the vehicle structure indicates 
that the instability is primarily due to the coupling with the Lnterpump mode. 
Therefore, this instability can be identified as an interpump mode type of 
instability. It was also found that the destabilizing influence on the thirty- 
fourth mode was also due to coupling with the interpump mode. 

3. 1'. 3. 2 Li ft off 

The results calculated for the liftoff condition are shown in figure 
11. In this case it is found that the fundamental mode (f ~ 2 Hz) is unstable 


If the structural gains given by the modal data had been employed, this mode 
would have remained stable with a minimum damping ratio of 0.003. 
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Figure 9d. Stability results for basic system: 
orbiter end-burn 
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Figure 9e. Stability results for basic system: 
orbiter end-burn 
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Figure 10, Unstable mode of basic system: 
orb iter end -burn 
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Figure 11c. Stability results for basic system 
liftoff 












if the first structural -mode frequency lies between about 2.2 to 2. 6 Hz; the 
damping ratio in this range reached -0.002. The other modes are found to 
be stable. The instability in the fundamental mode results from coupling 
between the first structural mode and the first feedline -type mode of the 
propulsion system (fig. 12). The inter pump mode primarily couples with 
the high-gain sixty-second structural mode (f 27 Hz). Although the coup- 
ling is significantly destabilizing, it does not produce an instability; the 
minimum calculated value of the damping ratio in this case being about 0. 002. 

3. 1. 3. 3 After SRB Separation 

The results calculated for this flight condition are shown in 
figure 13. In this instance the fundamental mode (f ~ 2 Hz) is seen to be 
unstable over a relatively wide range of the first structural mode frequency 
( —1. 9 Hz to ~2. 6 Hz). Again the instability results from coupling between 
the first feedline -type mode of the propulsion system and the first structural 
mode. The instability is severe, the damping ratio reaching about -0. 009. 
The high gain forty-sixth structural mode (f ~27 Hz) remains stable but loses 
a good deal of damping in the 22 to 26 Hz frequency range due to destabilizing 
coupling with the interpump mode; the minimum calculated value of the 
damping ratio in this instance being about 0. 003. 


3.2 System with Compliant Accumulator 


The effect on the system stability characteristics of the introduc- 
tion of a compliant accumulator are now examined. The employed accumu- 
lator has a 0. 05 7 (2 ft^) volume and is located at either the LPOP inlet, 

the LPOP discharge, or the HPOP inlet. 

3, 2. 1 Propuls ion -System Modes 

The introduction of the compliant accumulator into the propulsion 
system generally produced shifts of the feedline -type mode frequencies that 
were within the associated open-open and open-closed feedline frequency 
bands. The frequency of the interpump mode was not significantly changed 
by the presence of accumulators at either the LPOP inlet or at the LPOP 
discharge. However, the presence of an accumulator at the HPOP inlet 
reduced this frequency appreciably. For example, at the end-burn and liftoff 
conditions the reductions were from 24.4 to 14.9 Hz and from 23. 1 to 9.45 Hz, 
respe cti vely. 

3.2.2 Selection of Cases 


With the exception of the end-burn condition, the same set of 
structural modes was used as was employed previously in the analysis of the 
basic system. The exception at end-burn was the introduction of the twenty- 
first structural mode into the analysis for the accumulator at the HPOP-inlet 
location. This mode had a frequency of 15. 9 Hz and was introduced since it 
could possibly couple with the 14. 9 -Hz interpump mode. 
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3. 2. 3 


Stability Resalts 


3. 2. 3. 1 QrbLter End-Burn 

The results calculated for this case indicated that the stability of 
the lower structural modes used in the analyses (modes El through E3 0) was 
somewhat enhanced by the presence of the compliant accumulators at any of 
the three positions; the minimum value of the damping ratio calculated for 
these modes being in excess of 0. 0085. However, the case for the two 
proximate high-gain modes, the thirty-fourth and thirty-fifth modes, was 
radically different. The results for the thirty-fourth mode, shown in fig- 
ure 14a, reveal that the introduction of a compliant accumulator at the LPOP 
discharge renders the system unstable for structural mode frequencies in the 
range 22. 5 to 24. 6 Hz. The introduction of an accumulator at the LPOP inlet 
produces a decrease in damping with decrease in the structural mode fre- 
quency. The system does not become unstable, however, within the specified 
range of structural frequencies. In contrast to these results, an accumulator 
at the HPOP inlet is relatively beneficial. Similar trends are observed in the 
case of the thirty-fifth mode (fig* 14b). It is seen that the introduction of an 
accumulator at either the LPOP discharge or at the LPOP inlet intensified 
the instability that was previously predicted for the basic system. v The inten- 
sification is particularly severe for the LPOP discharge case where the 
minimum calculated damping ratio drops to about -0. 010 as compared to the 
worst case value of around -0. 001 for the basic system. In view of the poor 
showing of the LPOP-dis char ge accumulator at this end-burn condition, it 
was decided to drop this location from the subsequent analyses. Unlike the 
two LPOP locations, the introduction of an accumulator at the HPOP inlet is 
seen to be beneficial. The instability is eliminated and the minimum damping 
ratio maintained above 0. 006 over the specified range of structural frequen- 
cies. The pronounced effectiveness of the HPOP-inlet accumulator in the 
case of these two modes is undoubtedly due to the lowering of the interpump 
mode resonant frequency. The stability analyses indicated that the introduc- 
tion of possible coupling between the interpump mode (of the system with the 
HPOP-inlet accumulator) and the twenty-first structural mode was not sig- 
nificant - probably due to lower structural gain in this mode (about 2. 5 x 
10-5 1/kg). 


From practical consideration on the size of the accumulator, it is 
of some interest to see if the effectiveness of an accumulator at the HPOP 
inlet would be degraded by a reduction in the accumulator volume from the 
original 0. 057-rrP (2 ftp) value. This was checked bv performing calcula- 
tions for the thirty-fifth mode using a 0 , 029 -mP (l ftp) volume accumulator 
and a 0. 014-rrP (0.5 ftp) volume accumulator. The results calculated for 
these designs are shown in figure 15, where it is seen that the specified 
reductions in volume did not degrade the performance of the accumulator 
insofar as suppression of the interpump mode instability was concerned . 


Finally, it might be argued that the failure of the LPOP- inlet ac- 
cumulator to eliminate the instability in the high gain thirty- fifth structural 
mode was due to an inadequate level of compliance. To check this supposi- 
tion, an additional stability analyses was undertaken for this mode assuming a 

If the structural gain given by the modal data had been employed for the 
thirty- fifth mode, the system would have remained stable with the LPOP inlet 
accumulator and would have become unstable with the LPOP discharge ac- 
cumulator. In the former case the introduction of the accumulator was de- 
stabilizing in that the minimum damping ratio of the system was reduced. 
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Figure 14a. Stability results with compliant accumulator 
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Figure 15. Variation of stability with compliant accumulator 
volume for HPOP inlet location: orbiter end -burn 
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0. 57-m 3 (20 ft 3 ) volume accumulator at the LPOP inlet. This accumulator 
volume change produced a ten-fold increase In the accumulator compliance. 
The resulting stability results were essentially identical to those obtained 
for the 0. 057-m 3 (2 ft 3 ) accumulator thereby indicating that the inadequacy 
of t he LPOP-inlet suppressor lay with its location rather than in its size . 

3. 2. 3. 2 Liftoff 


The results calculated for the liftoff condition indicated that the 
introduction of an accumulator at either the LPOP inlet or the HPOP inlet 
would eliminate the instability previously predicted in the fundamental mode 
(fig. 16). The effectiveness of the LPOP-inlet accumulator in this case is 
attributed to the fact that the fundamental mode instability at liftoff involves 
coupling between a structural mode and the fundamental feedline -type mode 
of the propulsion system (unlike the end-burn instability that involves the 
Lnterpump mode of the propulsion system). It is also interesting to note that 
the HPOP-inlet accumulator, which was effective in eliminating the inter - 
pump mode type of instability at end-burn, is also effective in eliminating the 
feedline -type instability at liftoff. As was done at the end-burn condition, 
the change in the effectiveness of the accumulator at the HPOP inlet with 
decrease of the accumulator volume to 0. 029 m3 (1 ft3) and 0. 014 m 3 (0. 5 ft ^ ) 
was checked. The results in this case (fig. 17) were quite different and indi- 
cated a sensitivity to accumulator volume over the variation that was con- 
sidered. In the case of the 0. 029 -m 3 (1 f t ^ ) accumulator, the minimum value 

057 m 3 
mini - 

mum damping value became negative at the lower end of the structural fre- 
quency range. Thus, an accumulator of this size did not eliminate the feed- 
line -type mode instability. This failure is attributed to an insufficient lower- 
ing of the frequency of the fundamental propulsion-system mode. 

3 3 

With regard to the effect of the 0. 057-m (2 ft ) accumulator on 

the higher structural modes, some destabilizing influences on the fifty-first 
mode (f ~ 22 Hz) were found with the introduction of an accumulator at either 
the LPOP inlet or the HPOP inlet (fig. 18a). The mode remained stable in 
both cases, however, maintaining damping ratios in excess of 0. 005. A 
destabilizing influence of the LPOP inlet accumulator was exhibited with the 
sixty-second mode (f 27 Hz); the system barely remaining stable at the 
lower end of the frequency range of interest (fig. 18b). In contrast, the 
HPOP inlet accumulator had a beneficial effect on this mode increasing the 
minimum damping ratio from around 0. 002 for the basic system up to about 
0.006. 

3. 2. 3. 3 After SRB Separation 

As in the liftoff case, it was found that an accumulator at either the 
LPOP inlet or HPOP inlet will eliminate the fundamental mode instability that 
had been predicted for the basic system (fig. 19). This was particularly 
encouraging since the degree of instability in the fundamental mode for this 
flight condition was considerably more severe than for the liftoff event. A 


of the damping ratio was reduced to 0. 009 (from about 0. 013 for the 0. 
accumulator). In the case of the 0.014-m 3 (0. 5 ft 3 ) accumulator, the 
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Figure 16. Fundamental mode stability with compliant 
accumulator: liftoff 
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check was made on the sensitivity of the HPOP-Lnlet accumulator to reduc- 
tions in the accumulator volume. The results for a 0. 02 9 -m3 (1 ft^ ) ac- 
cumulator and a 0. 014-m^ (0. 5 ft^) accumulator are shown in figure 20. It is 
seen that the sensitivity is less than at the liftoff condition. There are re- 
ductions in the damping, but the system remains stable with damping ratios 
in excess of 0. 007. With regard to the high-gain forty-sixth mode (f ~ 27 Hz), 
the results (fig. 21) indicated that use of a 0, 057-m^ (2 ft^) accumulator at 
the HPOP inlet was beneficial and that use of the LPOP-inlet accumulator 
produced a decrease in damping with decrease in the structural mode fre- 
quency. These trends are similar to those observed in the high-gain struc- 
tural modes for the other two flight events. 


3.3 System With Resistive Accumulator 

3 3 

The design requirements on the 0. 057-m {2 ft ) volume resistive 

accumulator were Real Y_ ^ 0. 06 m^/MN s (1 Ln.^/sec) for 2 5 f < 30 Hz. 

a 

The specific values of the accumulator parameters that were selected to 
meet this requirement were 


R =1.7 MN s/m^ (0. 1 sec/in. 2 ) 

St 

L = 0.023 MN s 2 /m 5 (0. 00133 sec 2 /in, 2 ) 

As noted previously, the selected values were judged to be practically 
achievable on the basis of Delta Stage I vehicle experience (ref. 1). In addi- 
tion to this particular design, three other resistive accumulator designs 
were briefly examined to determine the sensitivity of the stability results to 
change in the accumulator parameters. 

3.3.1 Propuls Lon -System Modes 

The propulsion-system modes for the various resistive accumula- 
tor configurations are quite similar to the modes that were previously cal- 
culated for the corresponding compliant accumulator cases, the main dif- 
ferences in the results being an increase in the damping of the feedllne -type 
modes for the case of the LPOP-inlet location. 

3. 3. 2 Selection of Cases 


With the exception of some of the less interesting modal cases at 
liftoff, the same set of stability cases was employed as was used in the study 
of the compliant accumulator. The deletion of some of the modes at the lift- 
off condition (the eighteenth, twenty -sixth, and forty -ninth modes were 
dropped) was done to reduce the proliferating number of stability cases. 
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Figure 20, Variation of stability with compliant accumulator 
volume for HPOP inlet location: after SRB 
separation 









3. 3. 3 


Stability Results 


3.3,3. 1 Orbiter End-Burn 

In the case of the first through twenty -first structural modes, the 
results obtained with the selected resistive accumulator design were found to 
be essentially identical to the corresponding results obtained with the com- 
pliant accumulator. In the case of the thirtieth (f — 22 Hz) structural mode, 
there were some slight differences between the resistive and compliant 
accumulator results, the differences being more pronounced for the LPOP- 
inlet location. 

The results for the thirty-fourth and thirty-fifth modes (f '*-27 Hz) 
with a resistive accumulator at the LPOP inlet are shown in figures 22a and 
22b. Comparison of these results with the corresponding curves for a com- 
pliant accumulator at the LPOP inlet (figs. 14a and 14b) show a pronounced 
difference, in that the variation of damping ratio with structural frequency 
is no longer monotonic. Since the stabLlity picture in these modes is largely 
controlled by the interaction between the structural mode and the interpump 
mode of the propulsion system, it is believed that the differences in the re- 
sistive and compliant accumulator results are due to changes on the effective 
boundary condition on this interpump mode at the LPOP. In contrast to the 
LPOP- inlet location, the results obtained for a resistive accumulator at the 
HPOP inlet (figs. 22c and 22d) are very similar to those obtained with the 
compliant accumulator (figs. 14a and 14b). 

The net result of these resistive accumulator analyses was to indi- 
cate that the stability picture was essentially unchanged from that predicted 
for the compliant case . The lower structural modes remained stable; the 
instability in the thirty-fifth mode was eliminated by use of the accumulator 
at the HPOP inlet but was not eliminated by use of a resistive accumulator 
at the LPOP inlet. 

In this latter case, the suggestion could be raised that the ineffec- 
tiveness of the resistive accumulator at the LPOP inlet was peculiar to the 
particular design used in the analysis. To investigate this possibility, 
stability analyses were undertaken with the following additional resistive 
accumulator designs. 


Real Y 

a 


2 0. 06 m /MN s 


L = 0. 023 MN s 2 /m 5 

d 

R = 15 MN s/m 5 

3 


L = 0. 034 MN s 2 /m 5 

3 

R =8.5 MN s /m 5 
a 


(b) 
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Figure 22c, Stability results with resistive accumulator: 
or biter end -burn 
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Real Y > 0, 3 m 5 /MN s 
a 


L = 0. 007 MN s 2 /m 5 
a 

R =1.7 MN s /m 5 
a 


The results obtained with these accumulators at the LPOP inlet were similar 
to the previous results for the thirty -fifth mode and the instability was not 
eliminated. Thus, the results indicated that the ineffectiveness of the ac- 
cumulator at the LPOP inlet was due to the location rather than to the par- 
ticular accumulator design. 

Finally, a limited study was made to see if reduction in the volume 
of the resistive accumulator, located at the HPOP inlet, would reduce its 
effectiveness in eliminating the thirty-fifth mode instability. Calculations 
were made for a 0. 028-m-^ (1 ft^) volume accumulator and a 0. 014-m^ 

(0. 5 ft 3 ) accumulator. The resulting damping ratios were very close to the 
results for the 0. 057-m^ (2 ft3 ) volume accumulator thereby indicating that 
these volume reductions did not degrade the effectiveness of the accumulator 
insofar as elimination of the interpump mode instability was concerned. 

3. 3. 3. 2 Liftoff 


The results for the liftoff condition showed the same general 
trends as at the orbiter end-burn condition. The calculated damping ratios 
for the fundamental modes were close to the corresponding compliant ac- 
cumulator results. In the higher modes, the fifty-first (f ~ 22 Hz) and sixty- 
second (f ~ 27 Hz), pronounced differences were found between the resistive 
and compliant cases for an accumulator at the LPOP inlet. The results in 
the higher modes for the HPOP inlet location were similar to the correspond- 
ing compliant accumulator results. For this location, the effect of reducing 
the accumulator volume to 0. 028 m3 (1 ft^j and 0. 014 (0. 5 ft ^ ) was also 

checked. The results were almost the same as the corresponding curves 
obtained for the compliant accumulator at the HPOP inlet (fig. 17) indicating 
that the effectiveness of this device in eliminating the feedline -type insta- 
bility in the fundamental mode was significantly degraded by a volume reduc- 
tion to 0. 014 (0. 5 ft^). Again, this degradation in effectiveness is attri- 

buted to an insufficient lowering of the frequency of the first propulsion 
system mode. 

3. 3, 3. 3 After SRB Separation 

The same general trends were exhibited for this case. The results 
for the fundamental mode were very close to the compliant accumulator 
curves (fig. 19) while differences were found in the forty-sixth (f ~ 27 Hz) 
results for an accumulator at the LPOP inlet. 
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4. SUMMARY AND CONCLUSIONS 


This report provides a qualitative assessment of (1) the tendency 
for pogo instability of the Shuttle vehicle in the absence of suppression de- 
vices and (2) the effectiveness of accumulators in counteracting any tendency 
toward instability. Compliant and resistive accumulators are considered at 
three positions in the oxidizer system: LPOP inlet and discharge, and HPOP 

inlet. Stability analysis is performed by means of digitally computing eigen- 
solutions of the equations describing the coupled structure /propuls ion sys- 
tem. 


Engineering judgment was applied to construct a mathematical 
model for the study, with due consideration of the preliminary nature of the 
available data and the level of complexity necessary to deal with the primary 
aspects required for a stability assessment. Structural modes of vibration 
for coupled motion in the vehicle pitch plane were obtained from the Rockwell 
International Space Division. The propulsion system was simplified to the 
model shown in figure 3. The numerical values of the propulsion system 
parameters were based upon the primary geometric and dynamic performance 
features of the actual system as given by drawings and Rocketdyne's dynamic 
model of the main engine. The best available basis was employed to predict 
the cavitation compliances of the pumps, which is a method derived from 
empirical studies of measured compliances of the family of pumps used on 
Saturn vehicles. The analytical model is especially rigorous with regard to 
feedline hydraulics in that a continuous representation, including resistance 
effects, is employed to represent the multiplicity of damped organ -pipe 
modes. Contributions to the generalized force for a structural mode included 
frictional forces due to fluid resistance, fluid momentum-type forces, forces 
at turns in the fluid flow and at area changes, tank outflow effects, as well as 
engine thrust. 

In the course of the study some sixty-six basic combinations of 
flight condition /structural mode /accumulator were treated in order to develop 
comparative data. For the purpose of easy reference, the totality of these 
basic cases is presented in the following table together with an indication of 
whether the system was stable or unstable. In addition to these basic cases, 
the effects of particular variations in the compliant and resistive accumula- 
tor designs were also examined. 

The major findings from this study are as follows: 

1. The basic system can undergo instabilities involving a mode 
of the propulsion system which can basically involve either 
the in ter pump mode or a feedline mode. 

2. An instability involving the inter pump mode is suppressed 
effectively only by an accumulator located at the HPOP inlet. 
Such instability can be intensified by an accumulator at the 
LPOP inlet or discharge. 
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3. An instability involving a feedline mode can be counteracted 
effectively by accumulators at either the LPOP inlet or the 
HPOP inlet. 

4. There is no significant difference in the effectiveness of the 

compliant and resistive accumulators at the HPOP-inlet 
location. Reduction of the volume of either type of accumu- 
lator at the HPOP inlet to below 0. 03 (1 ftp) significantly 

degraded its effectiveness in suppressing feedline -mode 
instabilities. This same reduction of volume did not signi- 
ficantly degrade the effectiveness of either type of accumula- 
tor in suppressing the interpump mode of instability. 



Suppressor Location 1 LPOP Inlet 


2 LPOP Discharge 

3 HPOP Inlet 
X - calculated instability 

NA - no analysis 
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"With due recognition of the preliminary nature of the analysis, 
the primary conclusion of this study is believed to be valid: effective pogo 
suppression can be accomplished using an accumulator at the HPOP-iniet 
location, but not at the LPOF locations. The determination of which of the 
two basic types of accumulators, compliant or resistive, is best must await 
more detailed studies. It is particularly significant that future studies deal 
with the physical separation of the accumulator from the HPOP inlet, as 
dictated by engine functional and geometric considerations. 


With respect to tentative design requirements, it appears that an 
accumulator volume of between 0. 03 m^ (1 ft*) to 0. 06 m^ (2 ft^) would be 
appropriate for the compliant device. The same volume range together with 
the condition Real Y ^0.06 m^/MN (1 in.^/sec) appears to be adequate for 
the resistive accumulator. 


5. FUTURE STUDIES 


It is recommended that the study of passive suppression devices 
within the engine be continued and that both compliant and resistive accumu- 
lators be treated. It is recommended that the analytical model be upgraded 
by 


1. Use of a three -engine representation of the propulsion sys- 
tem that incorporates the ducting between the engines 

2. Use of updated structural-mode data and representation of 
the gimballing degree of freedom of the engines 

3. Consideration of the physical separation between the accumu- 
lator and the pump inlet 

Two additional factors that will require consideration at a later stage of the 
pogo stability analysis are the three-dimensional character of the interpump 
ducting and the effect of mechanical resonances of this ducting on the stability 
picture. Recently, it has been tentatively concluded that mechanical reso- 
nances on similar pump discharge ducting of the Delta Stage I vehicle are 
contributing significantly to observed instabilities. 

In the event that a physical design of a resistive accumulator is 
undertaken, analytical and experimental studies will be required to ensure 
linear resistive behavior; nonlinear resistive behavior was encountered with 
the orifice d liquid flow in the accumulator employed on the Delta Stage I 
vehicle. An alternative is the nonlinear resistive device which would entail 
accumulator testing within the operating engine environment to measure the 
random flow fluctuations into the accumulator and their effect on the accumu- 
lator resistance. This random flow has a considerable bearing on orifice 
resistance for small superimposed sinusoidal flow fluctuations. 
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APPENDIX A 


SYSTEM EQUATIONS 


A. 1 Fluid Dynamic Equations for Propulsion System 


The fluid dynamic equations for the individual elements of the propul 
sion system model are given below. 

First Feedline Segment 

P 1 = *11 P t + a iZ Q t + Q 'l3 X il 
Q 1 = a 21 P t + *22 Q t + °23 X il 
First Feedline Corner 


P = P 
1 2 


°2 A 2 Z 1 = Q 1 + A l x l 


Second Feedline Segment 


P 3 a ll P 2 + °T2 Q 2 " ^13 Z i2 
Q 3 = a 2l P Z + 3 22 Q 2 ' S 23 Z i2 


Second Feedline Corner 


P = P 
3 *4 


Q 4 + A 2 X 3 = °3 " A 2 Z 3 


Low-Pressure Pump 


Q 5 ' A 3 Z 4 = < Q 4 +A 2 X 4> + Q A1 - sC bl P 4 
P 5 = (tnpl)P 4 - Z pl Q 5 +R pl A 3 i 4 
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LPOP Discharge Suppressor 


Q 6 “ A 3 Z 4 " Q 5 " A 3 Z 4 + Q A2 

P 5 85 P 6 

Interpump Line 

P/ - P- = Z.Q, - A_R.z 
6 7 i 6 3 i 7 

Q, = Q 
w 6 7 


High- Pres s ure Pump 


Discharge Line 


Q 8 " A 4 Z 7 “ Q 7 " A 3 Z 7 + °A3 sC b2 P 7 


P 8 “ {rn 2 +l) P 7 “ Z p2 Q 8 + R p2 A 4 Z 7 


P 8 “ P 9 = Z d°8 " R d A 4 Z 7 


Q9 ~ Qg 


Injector and Chamber 


P 



Z.Q- - R . A z 
J 9 J 4 7 


P 

c 


= R c<°9 


V? 


) 


In the above equations, A^,Q. and P. denote the flow areas, absolute volume 
flows, and pressure perturbations at various locations within the system. 

The denote the volume flows from the pogo suppression devices: (m^+1) 
and (m 2 + l) are the LPOP and HPOP gains; and denote the cavitation 
compliance at the inlets to the LPOP and HPOP, respectively, while s is the 
Laplace variable; the x^, are the longitudinal and lateral velocities of the 
structure at various points of the system while and denote the aver- 
age translational velocities of the two feedline segments. It will be noted 
that the lateral velocity of the HPOP, thrust and injection chambers are 
taken to be the same. The individual impedances Z^ that appear in the above 
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equations are as follows: 


First Feedline Segment 


Z 


1 


LjS + R 


1 


Second Feedline Segment 


Z 2 = L 2 S + R 2 


Low Pressure Pump 



L , s 
P* 


+ R 


Pi 


High Pressure Pump 



L p2 s + R 


p2 


Interpump Line 

Z. = L.s + R. 

L l L 

Discharge Line 

z d = V + R d 

Injector 


Z. = L.s + R. 

1 J J 

The engine impedance that appears in the equations of motion for the com- 
bined str uctural/propuls ion system is defined as 


Z = L s + R 
e e e 
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where 


Li — L. -t L, + L — 

e j d p2 


R = R . + R + R . + R 0 
e j c d pZ 


The transmission coefficients for the first feedline segment are defined as 
follows : 


"n = "22 = cosh e i 


sinh 6 


*12 = " Z 1 8 , 


s inh 0 


O' = - A,R 


13 


11 0 , 


a 


21 


0 sinh 0 
Z i 1 1 


O' 


A i R i 


23 


1 - cosh 0 


where 6^ is defined by 


R 


e* = .22 1+ _y 

1 1 V sL 


with = ij/a where is the length of the feedline segment and a is the 
speed of sound in the liquid oxygen. The coefficients , a — , for the second 
feedline segment are s imilar with Aj , R j , Z ^ , 0^ be ing replaced by A^, 

Z^» respectively. 

The volume flow, Q^, from the pogo suppression devices are related 
to the propulsion system and specific suppressor parameters. The precise 
form of these relationships is given in Section A. 3. 


A. 2 Equation of Motion for Structural Mode 

The response of the vehicle structure is taken to be in the n^ normal 
mode. The structural motion x at some point ? is written 


st , 

x = q e 6 (r 
n n ' 
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where q n is the generalized coordinate and 9& n (r) is the associated mode 
shape. For the case that the structural system is defined with closed- 
bottom tanks, the equation of motion governing is 


M q + 21, <x> a + oj q 
n n n n n n m 


= ,^ n ( t b)6 R + (tb) + £F.. 


where ^(tb) is the modal tank-bottom pressure and Ls the relative volume 
outflow from the propellant tank. The detailed development of the contribu- 
tion of the outflow contribution is found in reference 3. From the results of 
reference 3, the tank-bottom pressure, Pj., is related to the vehicle motion 
by 


P 


t 


& 3 

n 



The Fj. comprises the drag forces on the feedline segments, the Lnterpump 
line and the discharge line, the forces at the two feedline corners, the forces 
on the LPOP and HPOP and the forces on the injector and thrust chamber. 
The precise form of these various forces is as follows: 

Drag on Feedline Segments 

-AiRi 

F i = ~ (p t “ p 2 + A i L i s x n ) 


(z) 


A 2 R 2 


(P, ~ p , 


" A 2 L 2 S Z f 2 


Drag on lnterpump Line 


F (z) 


A 3 R ~ A 3 Z 7^ 


Drag on Discharge Line 



A 4 R (Q 8 ‘ A 4 z 7 
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Forces at First Feedline Corner 


T7> ) — _ "PA 

F cl ' P 2 A 


f{ c1 = - P 2 A 2 - 1 


Forces at Second Feedline Corner 




F ( * } = p a a 7 + f 

c2 4 2 


= P.A + f 
c2 4 2 


(^) 


Forces at LPOP 


F p X l’ = ' P 4 A 2 ‘ f l 2 A^ + 


F pV = - P 5 A 3 




Forces at HPOP 


,(z) .. 

p2 


P 7 A 3 - P 8 A 4 + 2f 


(2i V 

l A 3' A 4 ( 


Force at Injector 


F. (z) = A.P q 
i 4 9 


Force at Thrust Chamber 


F t X) = A T P c 


In the above expressions f denotes the mean mass flow of the propellant 

7 = pVA 


71 



where p is the mass density of thejDropellant and V is the steady flow 
velocity. The items that involve f derive from the convective derivative 
in the fluid momentum equation. 


A. 3 Description of Accumulator 


The volume flow, Q^, from an accumulator in the system is written as 



Y P 
a 


where P is the pressure perturbation at the accumulator; Y a is the admittance 
of the accumulator. For a resistive device, it is convenient to express the 
admittance Y in the following form 

s/L 

Y = 5 

a 2 , 2 

S + ZL, w S t co 

a a a 


where L &f £ a and denote the inertance, damping ratio and natural fre- 
quency of the accumulator, respectively. In the special case of a purely 
compliant device, it is convenient to write 


where C is the compliance of the accumulator. 

E 


APPENDIX B 

SHUTTLE LOX SYSTEM PARAMETERS 


;tance, (1) 

MN s/m^ 

, 2 
(sec /in. 

R x = 

0.4 

(0. 0236) 

*2 = 

0. 11 

(0. 0064) 

R . = 
pi 

8. 8 

(0. 52) 

R 9 = 

p2 

44. 7 

(2.64) 


(1 


Resistance in SI 
in English units, 


units is based upon pressure divided by volume flow; 
weight flow is employed. 
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1.9 (0.11) 


R. = 

i 

R = 156 (9. 19) 

R c = 42.3 (2. 50) 

R d = 8.3 (0.49) 

R. = 60.2 (3. 56) 

INERTANCE, f2) MNs 2 /m 5 (sec 2 /in. 2 ) 
L } = 0.71 (0.0417) 

L 2 = 0.19 (0.0113) 



0. 017 

(0. 001) 

if 

ro 

0. 042 

(0. 0025) 

L. ~ 

L 

0. 22 

(0. 013) 

L 

e 

0. 38 

(0. 0225) 

L d = 

0. 25 

(0.015) 

L. 

0. 085 

(0. 005) 


CAVITATION COMPLIANCE 

See figure 4 

PUMP GAIN 

See figure 5 
TIMES (sec) 

Tj = 0. 0656 

Feedline Travel Time 

t 2 = 0.0178 


C bl 
C b2 « 
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AREAS (in. 2 ) 


A 1 = 

75. 7 

a 3 = 

31.2 

A 4 = 

12. 6 

A t = 

158 
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